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Abstract 

We report the first observation of the radiative hyperonic B decay PA7, using a 140 fb~^ 

data sample recorded on the T(4S) resonance with the Belle detector at the KEKB asymmetric 
energy e+e" collider. The measured branching fraction is B{B'^ — > pA'j) = (2.16^Q 53ib0.20) x 10~^. 
A search for pYP^ yields no significant signal, so we set a 90% confidence-level upper limit 

on the branching fraction of B{B^ — > pTP^y) < 3.3 x 10~^. 

PACS numbers: 13.40.Hq, 14.40.Nd, 14.20.Dh, 14.20.Jn 
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The 6 — i> S7 penguin diagram plays an essential role for the large rates of the observed 
radiative B —> K*^ l| decays. It is also a good probe of new physics beyond the Standard 
Model 0]. Recently, the Belle collaboration reported ^very stringent limit of O{10^^) on 
the branching fraction of two-body pA decays ^ but found an unexpectedly large 

rate for the three-body decay B^ pAvr" j^, which proceeds presumably via the b s 
penguin process. One interesting feature of the B^ pAvr" decay is that the observed 
proton-A invariant mass M^^ spectrum peaks near threshold. Naively, a suppression of 
0{aEM) is expected for the B~^ — > pA'j decay relative to B^ pA if the former process is 
bremsstrahlung-like. In contrast, the short-distance b ^ s'y contribution can lead naturally 
to a non-bremsstrahlung-like energetic photon spectrum and an enhancement of M^j^ at low 
mass; the former distribution can be compared to the recently measured b —>■ s'j inclusive 
photou energy spectrum g These features mot.vate our study of ^ pAr Some 
theoretical predictions [y] for the branching fraction of B^ pA'~) are at the 10 ^ — 10 ^ 
level, which is in the sensitivity range of the B-factories. 

We use a data sample of 152 x 10^ BB pairs, corresponding to an integrated luminosity of 
140 fb~^, collected by the Belle detector at the KEKB [tI asymmetric energy e^e~ collider. 
The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three-layer 
silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-fiight scintillation 
counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect mesons and to identify 

n 

muons. The detector is described in detail elsewhere j8|. 

To identify the charged tracks, the proton (Lp), pion (L^r) and kaon {Lk) likelihoods are 
determined from the information obtained by the tracking system (SVD+CDC) and the 
hadron identification system (CDC+ACC+TOF). Prompt proton candidates must satisfy 
the requirements of Lp/ {Lp + Lk) > 0.6 and Lp/{Lp + Lt,) > 0.6, and not be associated with 
the decay of a A baryon. The prompt proton candidates are also required to satisfy track 
quality criteria based on the track impact parameters relative to the interaction point (IP). 
The deviations from the IP position are required to be within 0.3 cm in the transverse (x-y) 
plane, and within ±3 cm in the z direction, where the z axis is opposite the direction of the 
positron beam. Candidate A baryons are reconstructed from two oppositely charged tracks. 
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one treated as a proton and the other as a pion, and must have an invariant mass within 
5(7 of the nominal A mass, as well as a displaced vertex and flight direction consistent with 
a A originating from the interaction point. To reduce background, a Lp/{Lp + L^^) > 0.6 
requirement is applied to the proton-like track. Photon candidates are selected from the 
neutral clusters within the barrel ECL (with polar angle between 33° and 128°) having 
energy greater than 500 MeV. We discard any photon candidate if the invariant mass, in 
combination with any other photon above 30 (200) MeV, is within ±18 (±32) MeV/c^ of the 
nominal mass of the 7r° (77) meson. The above selection criteria are optimized using Monte 
Carlo (MC) simulated event samples. 

□ 

Candidate B mesons are formed by combining a proton with a A and a photon j9||, 
each defined using the above criteria, and requiring the beam-energy constrained mass, 
Mbc = ^-Ebeam ~ Pb^ energy difference, AE = Eb — -Ebcam, to lie in the ranges 5.2 

GeV/c2 < Mbc < 5.29 GeV/c^ and -0.2 GeV < AE < 0.5 GeV. Here, pb and Eb refer to 
the momentum and energy, respectively, of the reconstructed B meson, and -Ebcam refers to 
the beam energy, all in the T(4S) rest frame. 

The dominant background is from continuum e^e~ — * qq processes (where q = u,d, s, c) 
while the background from B decays is negligible except for possible feed-down events from 
— > pS°7. This is confirmed with an off-resonance data set (10 fb~^) accumulated at an 
energy that is 60 MeV below the T(4S), and an MC sample of 120 million continuum events. 
In the T(4S) rest frame, continuum events are jet-like while BB events are spherical. We 
follow the scheme defined in Ref. Il0| and combine seven shape variables to form a Fisher 
discriminant |11|] in order to maximize the distinction between continuum processes and 
signal. The variables used have almost no correlation with Mbc and AE. Probability 
density functions (PDFs) for the Fisher discriminant and the cosine of the angle between 
the B flight direction and the beam direction in the T(4S) frame are combined to form the 
signal (background) likehhood Cs (-^b)- We require the likelihood ratio 71 = Cs/{Cs + C\,) 
to be greater than 0.75; this suppresses about 86% of the background while retaining 78% 
of the signal. The selection point is determined by optimizing Ns/ y/Ng + Nb, where Ng and 
Nfy denote the number of signal and background; here a signal branching fraction of 4 x 10~^ 
is assumed. 

We perform an unbinned extended maximum likelihood fit to the events with —0.2 GeV< 
AE < 0.5 GeV and Mbc > 5.2 GeV/c^ to determine the signal yields. The extended 
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FIG. 1: The distributions of AS (for Mbc > 5.27GeV/c2) and Mbc (for -0.135 GeV < AE < 
0.074 GeV) for pK-f candidates having M^^y < 2.4GeV/c^. The sohd, dotted and dashed 

hnes represent the combined fit result, fitted background and fitted signal, respectively. 



likelihood function is defined as 



-(AfA+TVs+Af,,) N 



i=l 



where Pa, Pe, and Pqq are the PDFs for pAj, pS°7, and continuum background, respectively, 
and A^A, A^s, and Ngg are the corresponding number of candidates. 

The pA'-f and pS°7 PDFs are two-dimensional smooth histograms determined by MC 



simulation. We use the parametrization first suggested by the ARGUS collaboration 



f{M^,) oc MbeV^l - (Mbc/^beam)'exp[-e(l " (Mbc/^bcam)')] , to model the background Mbc 
distribution and a quadratic polynomial for the background AE shape. We perform a two- 
dimensional unbinned fit to the AE vs Mbc distribution, with the signal and background 
normalizations as well as the continuum background shape parameters allowed to float. 

The AE distribution (with Mbc > 5.27 GeV/c^) and the Mbc distribution (with -0.135 
GeV< AE < 0.074 GeV) for the region M^^^ < 2.4 GeV/c^ are shown in Fig. [T] along with 
the projections of the fit. The two-dimensional unbinned fit gives a pA'y signal yield 

of 34.llg;g with a statistical significance of 8.6 standard deviations and a P+ pTP'y yield 



of 0.0 ±4.7. The significance is defined as y — 21n(Lo/Lmax), where Lq and Lmax are the 
likelihood values returned by the fit with signal yield fixed at zero and its best fit value, 
respectively {l3 |. 

We measure the differential branching fraction of pA'-j by fitting the yield in bins of M^^, 
as shown in Fig. |2l and correcting for the corresponding detection efficiency as determined 
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from a large MC sample of events distributed uniformly in phase space. The results of the fits 
along with the efficiencies and the partial branching fractions are given in Table m (In these 
fits, the signal yields are constrained to be non-negative.) The observed mass distribution 



in Fig. 

decays 



peaks at low pA mass, a feature seen also in —>■ pAn and ppK^ 



li. 



The photon energy spectrum in the T(4S) rest frame is measured using the same fit pro- 
cedure and is shown in Fig. 01 The predicted shape for the generic 6 ^ 57 process, obtained 
by MC simulation, is shown in Fig. El as the shaded histogram. The two distributions agree 
with each other within the statistical uncertainty of the present measurement, although the 
mean of the pA'j spectrum is slightly higher. 




3.5 

(GeV/c^) 



FIG. 2: The differential yield for B° pA-f function of Mp^. The shaded distribution is from 
a phase-space MC simulation with area scaled to the observed signal yield. 

We also study the proton angular distribution of the baryon pair system for Mp^ < 4.0 
GeV/c^. The angle 6x is measured between the proton direction and the 7 direction in 
the baryon pair rest frame. Fig. |31 shows the efficiency corrected B yield in bins of cos^x- 
This distribution confirms the b —>■ sj fragmentation picture where the A tends to emerge 
opposite the direction of the photon. 

The systematic uncertainty in particle selection is studied using high statistics control 
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TABLE I: The event yield, efficiency, and branching fraction (B) for each bin. 



M T (GeV/c2) Signal Yield Efficiency(%) B (IQ-^) 



< 2.2 


22.7t|i 


10.6 


1 41+0.40 


2.2 - 2.4 




9.8 




2.4 - 2.6 




9.3 




2.6 - 2.8 




9.9 


00+° °^ 
'-'•'-'U_o.06 


2.8 - 3.4 


0+3-^ 

U.U_3 4 


9.6 


U.UU_Q 23 


3.4 - 4.0 




9.6 





Total 33.81^;',' - 2.16+; 
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FIG. 3: Signal yield fraction versus photon energy in the T(4S) rest frame. The shaded histogram 
is the prediction from the inclusive 6 — > S7 MC simulation. 

samples. Proton identification is studied with a A ^ pn^ sample. The tracking efficiency is 
studied with a D* sample, using both full and partial reconstruction. Based on these studies, 
we sum the correlated errors linearly and assign a 4.1% error for proton identification and 
4.9% for the tracking efficiency. 

For A reconstruction, we have an additional uncertainty of 2.5% on the efficiency for 
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FIG. 4: Differential branching fraction versus cos9x in the baryon pair system. 

off-IP track reconstruction, determined from the difference of A proper time distributions 
for data and MC simulation. There is also a 1.2% error associated with the A mass selection 
and a 0.5% error for the A vertex selection [3]. Summing the errors for A reconstruction, we 
obtain a systematic error of 2.8%. 

The 2.2% uncertainty for the photon detection is determined from radiative Bhabha 
events. For the vr" and i] vetoes, we compare the fit results with and without the vetoes; the 
difference of the branching fraction amounts to 0.5%, and we quote this as the associated 
systematic error. 

Continuum suppression is studied by changing the selection criteria on TZ in the interval 
- 0.9 to see if there is any systematic trend in the signal fit yield. We quote a 2.5% error 
for this. 

The systematic uncertainty from fitting is 2.2%, which is determined by varying the 
parameters of the signal and background PDFs by ±lcr. The MC statistical uncertainty 
and modeling with six M^^ bins contributes a 4.4% error (obtained by changing the M^^ 
bin size). The error on the number of total BB pairs is 0.7%. The error from the sub-decay 
branching fraction of A — p7i~ is 0.8% jl^l ]. 

We combine the above uncorrelated errors in quadrature. The total systematic error is 
9.2%. 
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We see no evidence for the decay pTP-^. We use the fit results to estimate the 

expected background, and compare this with the observed number of events in the pTP"^ 



signal region (-0.2C 
limit on the yield 



GeV< < 0.04 GeV and Mbc > 5.27 GeV/c^) in order to set an upper 
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16 



17|. The estimated background for M^^ < 2.4 GeV/c^ (Mp^ < 4 



GeV/c^) is 45.5±6.3(84±0.2), the number of observed events is 44 (96), and the systematic 
uncertainty is 9.2%; from these, the upper limit yield is 14.9 (35.5) at 90% confidence level. 
The efficiency, estimated from the phase space MC simulation, is 12.1% (7.0%), so the 90% 
confidence-level upper limits for the branching fractions are B{B^ pTP'-j) < 0.8 x 10^^ for 
MpA < 2.4GeV/c2 and B^B"^ ptP^i) < 3.3 x 10"^ for M^x < 4.0GeV/c2. 

In summary, we have performed a search for the radiative baryonic decays B^ P^li 
and pS°7 with 152 million BB events. A clear signal is seen in the pA7 mode, and we 
measure a branching fraction of B{B+ pA-f) = (2.16 tl];! (stat) ± 0.20 (syst)) x 10-^ 
which is consistent with the upper limit set by CLEO 3]- The yield of the B^ pTP'^ 
mode is not statistically significant, and we set the 90% confidence level upper limit of 
B{B^ pfP^) < 3.3 X 10-*^ for M^a < 4.0GeV/c2. 
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